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nology,  Pasadena,  California. 


I 


ABSTRACT 

\ 

By  employing  the  Lighthill  ^'ideal  dissociating  gas,r 
approximation,  the  influence  of  the  equilibrium  disso¬ 
ciation  of  a  diatomic  molecule  on  Brayton-cycle  per¬ 
formance  is  demonstrated.  For  low  temperature 
ratios  it  is  shown  that  the  use  of  a  suitably  selected 
molecule  results  in  a  significant  improvement  in  cycle 
thermal  efficiency. 
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INTRODUCTION 


In  an  earlier  paper  we  investigated  the  influence  of  molecular  vibration  on 
Brayton-cycle  performance  by  employing  the  harmonic  oscillator  approxima¬ 
tion  [1]  .  This  approximation  allowed  us  to  write  the  equations  describing  the 
cycle  performance  in  a  convenient  dimensionless  form,  thus  admitting  a  con¬ 
siderable  degree  of  generality  in  our  analysis.  In  this  paper  we  treat  the 
influence  of  equilibrium  gas  dissociation  on  Brayton-cycle  performance.  A 
convenient  dimensionless  form  is  achieved  by  employing  the  approximations 
introduced  by  Lighthill  [2],  An  objection  which  could  possibly  be  raised  in  our 
use  of  the  Lighthill  approximation  in  this  application  is  that  it  embodies  the 
assumption  that  the  contribution  of  molecular  vibration  to  the  specific  heat  is 
frozen  at  one-half  of  the  classical  value  for  complete  excitation.  In  our  earlier 
work,  however,  we  have  demonstrated  the  influence  of  molecular  vibration  on 
cycle  performance  to  be  insignificant. 

In  the  succeeding  sections  we  show  that  gas  dissociation  implies  the  exist- 

2 

ence  of  at  least  one  extremum  in  the  cycle  thermal  efficiency.  We  then  present 
the  thermodynamic  analysis  and  discuss  the  influence  of  gas  dissociation  on 
Brayton-cycle  performance. 


1  Numbers  in  brackets  designate  References  at  end  of  paper. 

2  This  extremum  should  not  be  confused  with  the  maximum  implied  by  the 
introduction  of  component  efficiencies,  i.e.  ,  turbine  and  compressor  efficiencies. 
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EXISTENCE  OF  EXTREMA 


Following  our  earlier  argument  [l]  it  may  be  shown  that  the  maxi¬ 
mum  thermal  efficiency,  ,  of  a  Brayton-cycle  for  wnich  the  specific  heat  is 

a  constant  is  given  by  an  expression  which  is  independent  of  the  specific  heat  and 
is  of  the  form 


=  f  K,  *L.  — 


Here,  r\^  and  are  efficiencies  of  the  compressor  and  turbine,  respectively; 
T  is  the  absolute  temperature;  and  the  subscripts  refer  to  points  on  the  cycle  as 
shown  in  Fig.  1.  Now,  for  a  dissociating  diatomic  gas  at  sufficiently  low  tem¬ 
perature,  the  degree  of  dissociation  is  negligible  and,  hence,  the  specific  heat 
is  constant  if  we  neglect  the  contribution  of  the  vibrational  mode  of  motion.  At 


sufficiently  high  temperatures  the 
gas  is  completely  dissociated  and, 
again,  the  specific  heat  is  constant. 
At  these  extremes  the  maximum 
thermal  efficiency  is  given  by  Eq.  (1). 
It  then  follows  that  at  least  one  ex¬ 
tremum  exists  for  intermediate  de¬ 
grees  of  dissociation.  One  is  thus 
moved  to  investigate  by  detail  anal¬ 
ysis  the  nature  of  the  extrema  with 
a  view  toward  possible  improve¬ 


ment  in  cycle  performance. 

Fig.  1.  Tempe rature- Entropy  (S) 

Diagram  for  the  Brayton-Cycle 
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THERMODYNAMIC  ANALYSIS 


Consider  a  diatomic  species,  A^  ,  that  dissociates  according  to  the  equation 

A2-  2 A  .  (2) 

If  a  represents  the  fraction  of  A^  dissociated,  then  the  Lighthill  “ideal  disso¬ 
ciating  gas"  approximation  [2]  gives 

a2  pd  -D/RT 

1 -  =  —  e  •  (3) 

1  -  a  P 


Here, 

3 

p^  =  a  constant  for  a  particular  gas, 
p  =  the  actual  gas  density, 

D  =  the  dissociation  energy  per  mole  of  the  particular  molecule  at  0°K, 

R  =  the  universal  gas  constant,  and 
T  =  the  absolute  temperature. 

It  is  now  convenient  to  introduce  ,  p ^  and  as  units,  of  tem¬ 

perature,  density,  pressure,  and  specific  internal  energy,  respectively,  where 
Tj  ,  Pj  ,  and  are  defined  by  the  relations 


(4) 


P 


d 


(5) 


(6) 


A  method  for  rapid  calculation  of 


is  given  in  the  appendix. 
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and  where  W  is  the  molecular  weight  of  A^.  In  these  units,  Eq.  (3)  becomes 


Q2  1  -1/T 

1 -  =  -  e 

1  -  a  p 


(7) 


and  the  perfect  gas  law  and  specific  enthalpy  (also  in  units  of  U^)  may  be  written, 
respectively,  as 

P  =  pT(l  +  a)  (8) 

and 

h  =  u  +  |>  =  (4  +  a)T  +  a  4  .  (9) 


Along  reversible  isentropic  paths  we  may  also  write 


dS 


^i.dP  =  0 

T  T 


(10) 


Equations  (8),  (9),  and  (10)  may  be  combined  to  give  the  following  differential 
equation  relating  the  temperature  and  pressure  along  isentropic  paths: 


P  dT 
T  dP 


1  +  o  +  - 
2 


4  +  a  +  1q(1 


P 

T 


(11) 


With  the  equations  above  we  can  now  determine  the  maximum  thermal 
efficiency  of  a  simple  Brayton-cycle  including  the  influence  of  gas  dissociation. 
Referring  to  Fig.  1,  the  thermal  efficiency,  »  i*  given  by 

- 7 - 

Implicit  in  this  statement  of  the  enthalpy  is  the  approximation  that  the 
vibrational  motion  is  frozen  at  one-half  the  maximum  classical  value  of  RT. 


\ 
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(12) 


nTH 


The  component  efficiencies  are  defined  as 


and 


(13) 


(14) 


For  a  given  initial  pressure,  Pj  ,  and  specified  values  of  T^/Tj  »  *1  and 

t|j.  ,  Eqs.  (7)  through  (9)  and  (11)  through  (14)  allow  us  to  determine  for  any- 

assumed  value  of  .  By  systematic  variation  of  T^/T^  ,  the  value  of 

T?/T ,  which  maximizes  the  thermal  efficiency  may  be  found.  We  denote  this 

maximum  of  the  thermal  efficiency  by  .  The  actual  calculations  are  best 

carried  out  numerically,  and  we  have  programmed  these  equations  for  the  IBM 

709  digital  computer.  The  nature  of  the  solutions  obtained  is  indicated  in  Figs.  2, 

3,  and  4,  where  for  all  cases  we  have  taken  ti  and  r|  to  be  0.85. 

c  t 


Fig.  2.  Variation  of  Optimum  Thermal 
Efficiency  with  Characteristic  Dissocia¬ 
tion  Parameter,  T.=  RT  .  .  /D 
1  sink 
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Fig.  4.  Variation  of  Pressure 
Ratio  with  Characteristic  Dis¬ 
sociation  Parameter 
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DISCUSSION 


In  Fig.  2  we  show  the  influence  of  gas  dissociation  in  cycle  performance 
in  terms  of  the  percentage  change  in  maximum  thermal  efficiency  as  compared 
to  the  maximum  thermal  efficiency  for  a  cycle  in  which  the  specific  heat  is  con¬ 
stant.  It  will  be  recalled  that  temperature  was  made  dimensionless  with  respect 
to  the  bond  dissociation  energy  and,  thus,  one  may  "move11  along  the  abscissa  in 
our  figures  either  by  changing  the  dissociation  energy  or  the  sink  temperature. 

We  shall  call  T^  the  characteristic  dissociation  parameter.  It  may  be  seen 
from  the  figure  that  a  considerable  improvement  in  cycle  efficiency  can  be 
achieved  by  a  judicious  selection  of  diatomic  molecules. 

The  influence  of  dissociation  is  seen  to  be  greatest  for  low  temperature 
ratios;  however,  low  temperature  ratios  are  not  devoid  of  practical  application, 
since  it  may  be  readily  shown  that  temperature  ratios  of  less  than  two  are  required 
by  those  space  power  plants  which  are  radiator  weight  limited. 

It  is  interesting  to  note  that  a  number  of  maxima  and  minima  are  exhibited 
by  our  solution.  The  reason  for  these  extrema  may  be  inferred  from  Fig.  3,  in 
which  we  display  for  one  of  the  cases  shown  in  Fig.  2  the  fraction  of  gas  disso¬ 
ciated  at  various  points  in  the  cycle.  The  numbers  on  the  curves  refer  to  the 
corresponding  points  on  the  cycle,  Fig.  1.  As  an  example,  the  first  maximum 
in  Fig.  2  may  be  explained  by  noting  that  the  gas  is  entering  the  turbine  (point  3) 
partially  dissociated.  For  a  small  range  on  the  abscissa  the  gas  leaves  the  tur¬ 
bine  completely  associated,  the  energy  of  association  appearing  as  turbine  work 
so  that  the  cycle  efficiency  increases.  Continuing  along  the  abscissa,  however, 
dissociated  gas  begins  to  appear  at  the  turbine  exit  (point  4).  The  remaining 
energy  of  association  carried  by  the  gas  is  lost  during  the  heat  rejection  part  of 
the  cycle,  and  the  efficiency  begins  to  fall. 

In  Fig.  4  is  shown  the  pressure  ratios  required  to  achieve  the  efficiencies 
shown  above.  These  rather  large  pressure  ratios  could  present  a  severe  re¬ 
striction  for  practical  application. 
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At  a  number  of  meetings  concerning  space  power  plants  that  the  authors 
have  attended  in  the  past,  a  frequent  suggestion  for  improving  performance  has 
been  the  use  of  equilibrium  dissociating  gases.  Indeed,  Prof.  Lighthill  makes 
this  suggestion  in  his  paper.  To  our  knowledge,  no  calculations  have  been  forth¬ 
coming.  In  this  paper  we  have  clearly  indicated,  for  diatomic  molecules,  the 
magnitude  of  improvement  that  can  be  expected,  and  it  appears  to  be  appreciable. 
Further  study  with  a  view  toward  application  is  suggested. 
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APPENDIX 


In  this  appendix  we  show  a  rapid  method  for  the  calculation  of  the 
parameter  introduced  by  Lighthill  [2]. 

For  the  chemical  reaction  given  by  Eq.  (2),  the  equilibrium  constant, 
Kp  ,  is  defined  as 


Kp 


(15) 


Here,  P^  and  P^  are  the  partial  pressures  of  atom  A  and  molecule  , 
respectively.  In  terms  of  a,  the  mole  fraction  of  A  is  2cl/(1  +  a),  and  of  A^ 
is  (1  -  q)/(1  +  a).  Hence,  Eq.  (15)  may  be  written 


Kp 


(16) 


where  P  is  the  total  gas  pressure. 

From  the  fundamental  thermodynamic  formula  [3], 

Kp  =  «-AF°/RT 


(17) 


AF°  being  the  change  in  free  energy  produced  by  the  reactants  when  they  are  in 
their  defined  standard  states. 

We  may  write,  of  course,  the  tautology 


„  -l/R  A[{F°  -  D)/T]  -D/RT 

l\p  -  e  e 


(18) 


where  D  is  the  bond  dissociation  energy,  as  before.  Equations  (16)  and  (18) 
along  with  the  perfect  gas  formula  then  combine  to  yield 
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1  -D/RT 
—  e 
P 


(19) 


2 

a 

1  -  a 


W  -1/R  A[(F°  -  D)/T] 

<  — -  e 

4RT 

j 


Apparently,  the  quantity  in  brackets  is  to  be  identified  with  p^  ,  i.e., 


W  -1/R  A[(F°-  D)/T] 
Pd  4RT  6 


(20) 


The  value  of  A[(F°  -  D)/T]  has  been  tabulated  as  a  function  of  temperature  for 
a  large  number  of  substances  [4].  As  Lighthill  has  pointed  out,  p^  shows  a  very 
weak  temperature  dependence  for  many  diatomics. 
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